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ALTERNATIVE METHODS TO MODEL FRICTIONAL 
CONTACT SURFACES USING NASTRAN 
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SUMMARY 

Elongated (slotted) holes have been used extensively for the integration 
of equipment into Spacelab racks. In the past, this type of interface has been 
modelled assuming that (i) there is no slippage between contact surfaces, or 
(ii) there is no load transfer in the direction of the slot. Since the contact sur- 
faces are bolted together, the contact friction provides a load path determined 
by the normal applied force (bolt preload) and the coefficient of friction. This 
paper examines three alternate methods that utilize spring elements, exter- 
nally applied couples, and stress dependent elements to model the contacted 
surfaces. Results of these methods are compared with results obtained from 
methods that use GAP elements and rigid elements. 

INTRODUCTION 

Elongated holes have been used in the design of Spacelab Experiment Equipment 
mounting provisions. This type of joint is employed where large tolerances are allowed 
in one direction of the hole for the ease of integration. A simple way to model these 
interfaces is to use RIGID elements with the assumption that two connecting grid points 
are not moving against each other when loads are applied. Another common method is to 
use RIGID elements with the degrees of freedom associated with the longitudinal direction 
of the elongated hole released. When using this method, it is assumed that the joint did 
not carry load in its slotted direction. 

Due to the assumption involved, neither method yields realistic results. GAP elements 
have been used to achieve better results. However, when large numbers of GAP elements 
are used in a complex model, an unreasonably large amount of computer processing time 
is required to solve the system and, hence, is not economical or practical. 

In this paper, three alternative methods are investigated. The first method employs 
a spring element with a spring rate equal to the maximum load at the connection divided 
by the gap length. By using a spring element at the connection, the nonlinear frictional 
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force developed at the gap is replaced by a spring force that increases proportionally with 
gap distance. The second method used an externally applied couple to represent the 
frictional forces at the two contacted surfaces. This method creates a local realization of 
frictional forces when two grid points, representing two contacted surfaces, move against 
each other. In the third method, two contacted surfaces are connected by elements with 
stress-dependent material properties. The piecewise linear static analysis rigid format is 
ultilized to solve element forces at these elements. 

A simple stowage container and four supporting columns were developed to represent 
a typical installation in a Spacelab rack. The container was integrated on supporting 
columns using the above methods. The models were run on a SUN workstation using 
CSA/NASTRAN. The results obtained from each method were then compared. 

THEORY AND BACKGROUND 

Figure la shows a typical elongated hole and Figure lb shows force vs. relative 
displacement of an elongated hole. The force vs. displacement curves of other elements 
are shown in Wini r ep 2 through 7. 

From inspection, the stress dependent material element is the more appropriate ele- 
ment to simulate elongated hole behavior because the GAP element, spring element and 
coupled element would generate some error. A rigid element with all three translational 
degrees of freedom coupled (R123) should be used only when the frictional force is greater 
than the force at the connection. A rigid element with degrees of freedom associated with 
the longitudinal direction of the hole released (R23) should be used when, the frictional 
force and the displacement in elongated direction are small. 
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ment of GAP element Fig 7 Force vs displacement of nonlinear element 
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ANALYSIS AND RESULTS 

A simple box is constructed using QUAD2 elements. It is supported by four columns 
which were modelled using BAR elements. A plot of the structure is shown in Figure 8. 
In the rear, the box is mounted on two rear columns at grid points 13, 24, 37, and 48 
using rigid elements with all three translational degrees of freedom coupled. In front, the 
box is mounted to the front left column at grid points 17 and 41 and to the front ng t 
column at grid points 20, and 44, using a different integrating method. Gravitational loads 
varying from 0.2g to 4.0g are applied to the model in the X-direction, which is par e o 

the slotted direction. 

The material and properties of QUAD2 and BAR elements are chosen such that the 
box is heavier and more flexible than the columns. Since the columns are much stiffer than 
the box, displacement of the grid point on the column is small as compared to displacement 
of the grid point on the box. Displacement of the rear panel is closed to the displacement 
of rear columns because the rear panel interface points are mounted to the rear columns 
by rigid elements with all three translational degrees of freedom coupled. However, at 
the from panel, the relative displacement of front panel and front columns are greatly 
dependent on the type of connecting elements used. It should also be noted that the single 
point constraint forces developed at the front columns are also dependent on the amount 
of load that has been transfered to the columns through the front connecting elements. 

Properties of each type of connecting element are chosen in such a way that its de- 
formation can follow the elongated hole displacement curve as close as possible. The gap 
length and frictional force of the contacting surface were chosen arbitrarily at 0.5 inches 

and 15.0 lbs respectively. 

Displacement of the left front interfaces (grid points 17 and 41 were identical and 
were tabulated in Table 1). The single point constraint force of the left front column (gnd 
points 1005 and 1006) are tabulated in Table 2. Table 3 shows the displacement of the 
right front interfaces (grid points 20, and 44) while Table 4 shows a single point constraint 
force of the right front column (grid points 1007 and 1008). 
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TABLE 1 DISPLACEMENT VECTOR OF GRID POINT 
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TABLE 3 DISPLACEMENT VECTOR OF GRID POINT 20 AND 44 

RIGHT FRONT COLUME 
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DISCUSION 


Examination of the displacement of the left and right front panels of the box and the 
load distribution over four supporting columns indicates that types R123 and R23 are the 
two extreme cases. For rigid R123 types, the load was distributed evenly over four columns; 
meanwhile, for rigid R23 types, the load was carried by the rear columns only, and the 
box acts as a cantilever beam. The exact load distribution pattern is somewhere between 
the boundary established by these two types of connections. With finer load increments, 
the displacement curve of a stress dependent material element can be made to match the 
displacement curve of an elongated hole and, therefore, used as a reference to evaluate 
performance of other types of connections* 

In general, static analysis with GAP elements and piecewise Unear analysis employ 
an iterative scheme for solution. -They are costly to run especially when models become 
complex. For the problem at hand, the GAP element did not yield a corrected solution 
when under tension load, meanwhile, stress dependent material elements did not perform 
well when under compression. Therefore, these two types of elements were not suitable for 
* iuodel with multiple loading cases where the directions of the load at the interfaces are 

often unknown. 

A couple type connection shows the same characteristics as the R23 type, with the 
displacement shifted due to externally appUed couple. It is also difficult to use when the 
direction of the loads at the connection is unknown. However, the solution can be obtaine 

with less computing time. 

The solution for rigid elements and spring elements always contains some degree of 
error, but can be obtained with less computing time, and can be used in problems involv- 
ing multiple loading conditions. The spring rate can be adjusted to control the relative 
displacement of the two connecting grid points. 

RECOMMENDATION 

II is recommended that for problem, with multiple loading conditions, model, should 
be run first using R123 type connections. Then, for interface, that develop forces larger 
than the frictional forces of the elongated hole, the type R123 connection should be re- 
placed with appropriate spring elements. This methodology yields reasonab e res s wi 
min imum computer run time. 
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